TPI-MINN-01-36 UMN-TH-2020-01 



Best Values for the CP-odd Meson-Nucleon Couplings from 

Supersymmetry 



Maxim Pospelov 



Theoretical Physics Institute, School of Physics and Astronomy 
University of Minnesota, 116 Church St., Minneapolis, MN 55455, USA 

and 

Centre for Theoretical Physics, CPES, University of Sussex, Brighton BN1 

9QJ, United Kingdom^ 



Abstract 

In the supersymmetric models, the dominant sources of the hadronic 
flavor-diagonal CP violation at low energy are the theta term and 
the chromoelectric dipole moments of quarks. Using QCD sum rules, 
we estimate the preferred range and the best values for the CP-odd 
meson-nucleon coupling constants induced by these operators. When 
the theta term is removed by the axion mechanism, the size of the 
most important isospin triplet pion-nucleon coupling is estimated to 
be g^NN = 2 x 10~ 12 (d u — dd), where chromoelectric dipole moments 
are given in units of 10 -26 cm. 
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The search for CP violation in flavor-conserving processes is of paramount 
scientific importance. The suppression of CP-violating effects induced by the 
complex phase of the Kobayashi-Maskawa matrix allows the use of electric 
dipole moments (EDMs) of neutrons or heavy atoms as well as T-odd asym- 
metries in the decays and scattering of baryons as powerful tools for probing 
new physics beyond the Standard Model. 

The wide separation between the energy scale of "new physics" (super- 
partners, technicolor, etc.) and the characteristic momenta of particles in 
non-accelerator experiments permits consideration of only the first few terms 
in the effective CP-odd Lagrangian. In the minimal supersymmetric models 
only the theta term, three-gluon operator, and EDMs and color EDMs of 
light quarks are important: 

5C= - V m q q(l + i6 ql5 )q + 6 G -^GG + wGGG 

q=u,d,s 

~2 Z d Q^ Fa ^q~^ Z d q qg s Gaj 5 q, (1) 

q=u,d,s q=u,d,s 

where GG = G%G%, Go = t a G%a, lu and GGG = f abc G^G b ua G^. The 
coefficients in flU) are generated by the CP violation in the SUSY breaking 
sector and evolved down to 1 GeV, which is the border line of viability of 
the perturbative quark-gluon description. 

In this Letter we present a systematic study of the transition from this 
Lagrangian to the effective T-odd meson-nucleon interactions which deter- 
mine the magnitude of the CP-violating nuclear moments and the T-odd 
asymmetries in nucleon scattering. T-odd nuclear forces are the main source 
for the EDMs of heavy diamagnetic atoms (see e.q. Jl|). The quality of 
constraints imposed on supersymmetric models from a recently improved 
measurement of the EDM of the xenon and mercury atoms 0, |3j (as well 
as of future experimental efforts with the EDMs of diamagnetic atoms and 
the T-odd nucleon scattering J4|) depends crucially on the treatment of QCD 
and nuclear effects, i.e. on the extraction of limits on di from the experimen- 
tal bound on the atomic EDM. The implications of this powerful constraint 
for the CP violation in the supersymmetric models have been emphasized 
in Refs. [||, ||], and numerically exploited in |7|]. The purpose of this work 
is to give "state-of-the-art" estimates for various T-odd nucleon-meson cou- 
pling constants, i.e to find their best values in terms of the coefficients in eq. 
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([I]). This problem is reminiscent of Ref. || which estimates various P-odd 
meson-nucleon couplings in terms of the Fermi constant. 

T-odd nuclear forces inside the nucleus can be approximated by a meson 
exchange with one of the meson-nucleon couplings being T- violating ||. It is 
natural to expect that pion exchange dominates in the T-odd channel. The 
coupling of nucleons with pions can be conveniently parametrized [K| |TTJ as 



£ CP = sXat^tW + £ ] NN NNn° (2) 

These couplings are generated by the theta term and by the color EDMs of 
quarks, g = g(9, di), where 9 = 9 G + J29 q . Couplings that change the isospin 
by two units can be generated only at the expense of an additional m u — 
suppression and are ignored in the present analysis. g^N^) * s rather well 



known [Hfl, as it can be deduced from the size of the (N\uu — dd\N) matrix 
element. For most of the models of CP violation including minimal SUSY 
models, 9 has to be removed by the Peccei-Quinn (PQ) symmetry leaving 
quark color EDMs as the dominant source for CP-odd nuclear forces. The 
contribution of the three-gluon operator GGG to g is additionally suppressed 
by m q and can be neglected. 

The first step in the calculation of g^ N (9, d u + d~d) and g^ N (d u — dd) 



is the reduction of the pion field by means of PCAC [13|, Fig. la. The 
smallness of the t-channel pion momentum compared to the characteristic 
hadronic scale justifies this procedure, 

(N7r a \0 C p\N>) = ^(N\[Ocp,JS 5 }\N')- (3) 

In 

The commutator of the zero component of the axial current with CP- violating 
operators Ocp = qg s Ga-y 5 q can be easily computed, leading to the matrix 
elements of the qg s Gaq operators over the nucleon state [|13|]. However, eq. 
(|3]) is an incomplete result. A second class of contributions was pointed out 
in Ref. and in Refs.[|14], in the context of the neutron EDM problem. 
They consist of the pion pole diagrams, Fig. lb, which contribute at the 
same order of chiral perturbation theory. Indeed, the quantum numbers of 
the qg s Gaq operators allow them to produce zero-momentum 7r°'s from the 
vacuum. The pion-nucleon scattering amplitude at vanishing pion momen- 
tum is proportional to the first power of the quark mass, whereas the pion 
propagator contains l/m q , so that Fig. lb and la both contribute at the 
same 0(d q m^) order. 



3 



(a) 



(b) 



x Qgp 

71 



N 



-X 



N 



Figure 1: Two classes of diagrams contributing to the CP-odd pion-nucleon 
coupling constant. 

Using the low energy theorems that relate the pion-nucleon scattering 
amplitude with the matrix elements oim q qq over the nucleon state, we arrive 
at the following intermediate result for the NNk vertex: 

——(N\d u (ug s Gau — rn^uu) — d d (dg s Gad — rn^dd)\N) 



oZ) i 2 | it | d | dg 

29 + m 1 1 



(N\uu-dd\N). (4) 



In this expression, m* = m u rrid/ {jn u + md) and m\ = (0\qg s Gaq\0) / (qq) = 
— (0.8 ± 0.1)GeV 2 [jnj parametrizes the strength of the quark-gluon dim=5 
vacuum condensate. In our case, this originates from the {7T \qg s Ga , j5q\0) 
matrix element and the minus sign is included into the definition of 
for convenience. An alternative way of obtaining the amplitude (f|) is to 
chirally rotate quark masses to the basis where pions cannot be produced 
from the vacuum, (7T°\2Yl,qTn q 9 q j5q + Yl,d^qg s Ga^q\Q) = 0, while keeping 
the theta term fixed, 9 q — const. This eliminates diagrams lb, but creates 
an additional contribution to la, leading to the same result (f|). 

When the PQ mechanism is activated, removing the theta term, the min- 
imum of the axion vacuum is shifted from 9 = by the color EDM operators 



16]. It turns out that the true minimum is such that the square bracket in 
eq. (B) is zero so that only the first line survives. This leads to a relatively 
simple expression for the couplings 

A = ^^W»-n*\p) (5) 

£L=^j^(N\H u + H d \N) 
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Figure 2: Various classes of diagrams contributing to the OPE for the f 
structure. 



in terms of matrix elements of the H u and Hd operators: 

H u = ug s Gau — rn^uu] = dg s Gad — rn^dd 



(6) 



Previously, using a combination of QCD sum rules and scaling arguments, 
Ref. |L3| estimated that (N\qg s Gaq\N) ~ ^ml(N\qq\N). Another analysis 



fj| finds similar estimate (N\qg s Gaq\N) ~ m^{N\qq\N) . Obviously, these 



?(°) 



and 9 ( nNN 



estimates are not sufficient to derive a reliable answer for g^ N 
because of the additional — m^(N\qq\N) contribution coming from diagrams 
lb. The danger of mutual cancelation between the two contributions was 
realized in Ref. where the need for a dedicated analysis of (N\H U ^\N) 
was emphasized. In the rest of this paper we derive the QCD sum rules 
|T9|j for the matrix elements of the H u u\ operators. The advantage of this 



approach is that the operator product expansion (OPE) will contain similar 
vacuum condensates for both sources, qg s Gaq and —vn^qq 
to trace possible cancelations. Following Refs 
generalized nucleon interpolating current, 



which allows us 
T5| , |20| , we introduce the 



(7) 



which combines the two Lorentz structures, j\ = 2e a b c {u r [C^db)u c and ji = 
2e abc (ulCd b )^u c . 

We compute the OPE for the correlator of this current in the presence of 
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d u (d)H u (d) external sources 

n(Q 2 ) = z / d 4 xe ip - x (0\T{ m (x)r] N (0)}\0) 



d u ,d d ' 



(8) 



where Q 2 = —p 2 , with p the current momentum. 

We limit our calculation to the Lorentz structure proportional to f be- 
cause it is less susceptible to direct instanton contributions and excited reso- 
nances than the chirally even structure proportional to 1. Relevant diagrams 
for this correlator are shown in Fig. 2. After a straightforward calculation, 
we find 

,2\ i 



* 2 U 



In 



'3v 

p 2 



7r LO + 4ln 
p z 



Mr, 



n NLO + 



:7T 



NNLO 



The leading order term is given by the diagrams 2a-2b, 



7T 



LO 



3m 2 , 
"64" 



d_(5 - 2{J - 3p 2 ) + d+(l - py 



(9) 



(10) 



Here we have introduced the combinations d + = d u + dd and d- = d u — da- It 
turns out that the diagrams 2a, where the external source enters through the 
(0\qpq\0) structure, give large and opposite sign contributions for the g s qGaq 
and —rri^qq sources so that their combined effect in H q is nil. Fortunately, 
this cancelation does not hold for diagrams 2b that give (|IUD. 
The next-to-leading term corresponds to diagrams 2c, 



7T 



NLO 



3d_ 
~32 



rog(l - (3 2 ) - g 2 (GG) 
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(7 - 2/3 - 5/3 2 ) 



d-i 



+^g 2 s (GG)(l-P) 2 , 



(11) 



which contribute to the OPE (Q) with the log of the infrared cutoff Air. 
g 2 {GG) ~ 0.4 — 1 GeV 4 is the vacuum gluon condensate. The next-to-next- 
to-leading order 

24 

"° (13 -2(3 -llf3 2 ) 



TV 



NNLO 



Xs n 2 (l + 2[3-3(3 2 ) 



3mk 



16 



d + 



XT*\i-py 



(12) 



+■ 



16 



-11 - 14/5 + (3 2 ) - 9 ^ GG "> (5p 2 + 8/3-1) 

24 
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contains the vacuum polarizabilities, 

XS) T=Jd 4 x(0\T{uu±dd{x),H u ±H d {0)}\0), (13) 

and the vacuum factorization assumption has been made in (|T2|). 

The sum rules prescription involves matching the OPE with the phe- 
nomenological part, n OPE (Q 2 ) = Tl phen (Q 2 ), where 



p 1 n phe 



^ 2 9-nNN m N 

[p 2 — m 2 N ) 2 ' p 



+ 



A 



m 2 N 



(14) 



contains double and single pole contributions, and the continuum. After 
Borel transformation of the sum rule n OPE (Q 2 ) = n phen (Q 2 ) we obtain 



(gg) 
M ' 



7T E 



71 



NLO 



M 2 



-(In 



M 2 
A 2 " 



-0.58) 



7T 



NNLO 



M 2 



exp 



' M 2 



(2\ 2 m N g nNN + AM 2 ) + M~ 2 B exp 



M 2 



(15) 



Here s is the continuum threshold and Eq = 1 — e _S//A/2 . A and B parametrize 
the contribution of excited states and are assumed to be independent of M. 

It is reasonable to start the numerical treatment from a simple estimate, 
a la Ioffe [21], which assumes the dominance of the ground state and the 



LO OPE term, and eliminates A using the nucleon mass sum rule for p 1 
Separating different isospin structures, we find 

3 4vr 2 |(gg)|m 2 
2 



iJ-wNN 



(d u - d d ) 



m N f w M 2 



Fx{P) 



-(o) 

9ttNN 



(du + d d )l^M^E M . 



(16) 
(17) 



10 m N f w M 2 

Here F x {p) = (5 - 2(3 - 3/? 2 )/(5 + 2(3 + 5(3 2 ) and F (J3) = 5(1 - (3) 2 /(5 + 
2(3 + 5/3 2 ). -Fi(O) = F (0) = 1. To get numerical estimates, we choose (3 = 0, 
extensively used in lattice simulations. It is well known that the ji current 
has a much better overlap with the nucleon ground state and Ai ^> A2. 
Substituting M = lGeV, we obtain 

_ 12 d u - d d I (qq) \ 



U-kNN 



-(0) 
9ttNN 



3 x 10" 



0.6 x 10" 



10" 26 cm (225MeV) 3 0.8GeV 2 
12 d u + d d \(qq)\ \m 2 Q \ 



10" 26 cm (225MeV) 3 0.8GeV z 



(18) 
(19) 
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In most SUSY models, d u u) = loop factor xM^ sy x a linear combination 
of m u and m^. When combined with (qq) from eqs. (p!8|-pj|), this forms m^f 2 
times a function of m M /m^, thus eliminating a major source of uncertainty 
in EDM calculations due to the poor knowledge of m u + [^4[ 15, 20f| . The 



estimate ( |l"8~D is twice smaller than the value of g^N use d in Q. Also in 
agreement with |J, eqs. ( |18HT9| ) suggest that g^NN/fJ^NN ~ 0.2d + /cL. 

For a more systematic analysis, one has to include NLO and NNLO 
terms in the OPE. Here, we immediately face the problem of the unknown 
vacuum condensates Xs,t- Even though the vacuum correlators (qq, qq) can 
be determined using chiral perturbation theory p2| , there is no direct in- 
formation on (qq, qg s Gaq) other than that it is likely to be comparable with 
m1(qq, qq). At this point we would like to take advantage of the possibility to 
choose (3 in such a way as to minimize higher order terms in the OPE. We note 
that xs m ( |12D is multiplied by 1+2/3— 3/3 2 which becomes at (3 = — 1/3 and 
1. The choice of (3 = 1 also suppresses the leading order, while (3 = —1/3 
maximizes it. For the expected size of xs @, Xs ~ ±0.16 x ml(qq) f~ 4 , 
we can choose (3 in a range such that the whole square bracket in front of 
cL in eq. ( |T2"D is zero. This gives a range of interpolating currents around 
£ = -1/3, 

- 0.5 < /? < 0, (20) 

where we can tune the NNLO terms to zero in the g^^ channel. Variation 
of j3 in this range contributes to an estimate of the uncertainty in our analysis. 
In the g^NN channel there is no obvious choice of (3 that would remove xt 
and leave the leading order term un-suppressed, so we will choose the same 
(3 as for g^NN- We note that this range is close to (3 = as used for ( P~S) - [TD| ) . 
One should also worry about the dependence on Aj^ in NLO. Remarkably, 
in the range (|2Ti|), this dependence is softened by cancelation of the tUq and 
g 2 s (GG) terms. 

The preferred range for g^ N and g^ N is determined according to the 
following procedure. We take the OPE side of ([TBj) at the lower point of 
the usual Borel window, M 2 = 0.8 GeV 2 , and vary it through the range 
of parameters -0.5 < (3 < 0, 300MeV < A IR < 500MeV, 0.7GeV 2 < 
|m 2 | < 0.9GeV 2 , 0.4GeV 4 < g 2 s (GG) < lGeV 4 , and 2GeV 2 < s < 3GeV 2 , 
0.8GeV 6 < (2tt) 4 A 2 < 0.9GeV 6 and finally -6 GeV" 1 < Xr/\(qq)\ < 6 
GeV -1 . On the r.h.s of ([15]) we assume the dominance of the double pole 
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contributions for M 2 = 0.8 GeV 2 and allow for a 50% correction due to the 
presence of the unknown parameters A and B, thus effectively widening the 
allowed range for g n NN- As expected, the couplings are most sensitive to the 
value of 771q. The final results are presented in Table 1. Our "best" value for 
g^NN is determined by averaging over (3 and choosing the central values for 
condensates, which also suppresses the logarithmic term. In order to sepa- 
rate the contribution of g^N from the A and B terms, we impose a relation 
among A, B and 7r i0 obtained by requiring the same large M 2 asymptotic 
behavior for both sides of (p~5|) . The resulting sum rule is fitted numerically 
and produces a result 1.5 times smaller than the naive estimate (fl8|) . For 



g^N N the best value cannot be determined as the OPE side changes sign 
depending on the value of xt- 

Also included in this table are the preferred ranges for the CP-odd cou- 
plings of nucleons with rj, p and u mesons. The couplings with p and u have 
the EDM-like structures —^N{d v V^ — d^V^a^^N with properly arranged 
isospin indices. They can be easily extracted from the calculation of the neu- 
tron EDM d n , induced by d u ^ d ) [Tj| after a simple re-assignment of charges 



for the external vector currents. Best values for g P NN and g^NN follow from 
the central values of d n (d u ,dd) given in fR)| . Finally, the coupling to the rj 
meson is dominated by the strange quark chromoelectric dipole moment d s 
in the isospin-singlet and by d u — dd in the isospin-triplet channels, and in 
both cases only the expected range can be quoted. 

In conclusion, we have shown that the size of the CP-odd pion-nucleon 
constant generated by quark chromoelectric dipole moments is given by the 
matrix element of qg s Gaq—m^qq over the nucleon state. We have constructed 
a QCD sum rule for this matrix element and determined the preferred range 
and the best value for the g^NN coupling. The upper part of the preferred 
range agrees with previous estimates. However, in the interpretation of the 
experimental limit on the EDM of the mercury atom |§ in terms of limits 
imposed on new CP-violating physics, a more conservative value g^N = 
2(d u — rfd)/10 _14 cm should be used. This translates the result of Ref. (see 
[0. for details) into the bound \d u — dd\ < 2 x 10 -26 cm. This constraint 
provide a sensitive probe of CP violation in the supersymmetric spectrum 
up to Msusy of few TeV. 

Acknowledgments The author would like to thank P. Herczeg, M. 
Voloshin and especially A. Ritz for helpful discussions. 
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coupling preferred range best value 

g^ N xlO 12 (1 to 6) d_, 26 2 cL 26 

^ xlO 12 (-0.5 to 1.5) d + , 2 6 

xlO 12 ~ (-1.5 to 1.5) 4,26 

9r]NN x 10 12 (-0.3 to 1) d_, 26 

g$ N (0.7 to 2) rf + 1.4 d + 

sJjvjv (0-5 to 1.5) 0.9 

sSvjv (-1-5 to -0-5) d+ -0.9 gL 

& (-2 to -0.7) d + -1.4 d_ 



Table 1: Preferred ranges and best values for the CP-odd meson- nucleon 
coupling constants induced by quark chromoelectric dipole moments. All 
values contain an overall multiplier |(gg)|/(225 MeV) 3 . gL,26, c?+,26 and (i S! 26 
are oL and <i+ and c? s in units of 10~ 26 cm. 
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